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ABSTRACT 

Photometric data from the literature is combined with triaxial mass models to de- 
rive variation in the intrinsic shapes of the light distribution of elliptical galaxies 
NGC720, 2768 and 3605. The inferred shape variation in given by a Bayesian prob- 
ability distribution, assuming a uniform prior. The likelihood of obtaining the data 
is calculated by using ensemble of triaxial models. We apply the method to infer the 
shape variation of a galaxy, using the ellipticities and the difference in the position 
angles at two suitably chosen points from the profiles of the photometric data. Best 
constrained shape parameters are found to be the short to long axial ratios at small 
and large radii, and the absolute values of the triaxiallity difference between these 
radii. 

The elliptical galaxies of our present investigation are very flat, with ellipticity 
typically around 0.3 or more. We find that the expectation values of the short to long 
axial ratio of these galaxies are around 0.5. 
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1 INTRODUCTION 

Intrinsic shapes of t he indiv i dual e ll iptical galaxi e s hav e 



been investigat ed by Binnev (Il985 h. Tenies et al 
IStatlerl lfl994aHbh . iBak and Statlerl (|2000h . IStatlei 



1993), 



2001), 



and Statler et al.l l|2004 ). These authors have used the kine- 
matical data and the photometric data, and have used the 
triaxial models with the density distribution p(m 2 ), where 
m 2 — x 2 + y 2 /p 2 + z 2 /q 2 with axial ratios p and q. Here, 
(x,y,z) are the usual Cartesian co-ordinates, oriented such 
that x-axis (z-axis) lies along the longest (the shortest) axis 
of the model. It was shown analytically that the projected 
density of such a distribution p(m 2 ) with co nstant (p, q) 
is stratified o n simila r and c o-aligned ellipses ([Stark 1977; 
lBinnevlll985l ). IStatlerl l|l994al ) uses (apart from the kinemat- 
ical data) a constant value of ellipticity, which is an average 
over a suitably chosen range of radial distance, for the shape 
estimates. The shape estimates are robust, and are described 
by a pair of the shape parameters, namely the short to long 
axial ratio c\l of the light distribution and the triaxiality Tm 
of the mass distribution. 

A complementary problem was attempted by 
ijChakrabortv et al.ll2008l . hereafter C08), wherein variation 
in the intrinsic shapes of the light distribution of ellip- 
tical galaxies was investigated by using triaxial models, 
which exhibit ellipticity variation and position angle twist. 
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These models are fixed by assigning the values of axial 
ratios (po,qo) and (poo,goo) at small and at large radii, 
respectively. These axial ratios are related to triaxialities 
To and Too, respectively, at small and large radii. We use 
Bayesian statistics, and obtain the variati on in the shape , 
following the methodology described in IStatlerl (|l994al ). 
We find that the marginal posterior density (MPD) is 
likelihood dominated, so that it relatively insensitive to the 
unknown prior density. We use a flat prior. We use a large 
ensemble of models, so that the shape estimates may be 
model independent. 

The basic ingredients of our method are the same as in 
IStatlerl (|l994al ), and we adopt all the necessary alterations 
described in C08. We use (qo, To, goo, Too) as the shape pa- 
rameters and use the ellipticities e;„ , e out and the position 
angle difference O ou t — Bin at two suitably chosen points 
Ri n and R ou t from the profiles of the photometric data of 
the galaxies. We find that the best constrained shape pa- 
rameters are qo , qoo and the absolute value of the triaxiality 
difference T d , defined as \T d \ = \T X - T |. 

C08 have estimated the shapes of 10 elliptical galaxies 
which are comparatively rounder, with ellipticities ^ 0.3. 
We now investigate shapes of three more galaxies, namely 
iVGC720, 2768 and 3605. These are very flat galaxies with 
ellipticity around 0.3 or more. We find that the expectation 
values of the short to long axial ratio of these galaxies are 
around 0.5. We use triaxial models which are very flat. We 
take models with the lower limit of (go,9oo) ~ 0.3 for our 
shape investigation. We find that a class of very flat triaxial 
models develops several undesirable features (sect. 2 and 
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Appendix A), and are not employed in the present shape 
estimates. 

Determination of the intrinsic shape using photometry 
is important because the number of galaxies with good pho- 
tometric is many more than those with good kinematics. 
Besides, the results obtained by alternative models and tech- 
niques can be used for a comparison. Photometry constrains 
the flattening (go, ?oo) but can not constrain (To, Too). Thus, 
our work is complementary and not contradictory to that of 
Statler and his coworkers. 

Sect. 2 presents the models. The necessity for the choice 
of small values of the lower limits of (go, goo), and the intrin- 
sic shapes of the galaxies are presented in sect. 3. Sect. 4 is 
devoted to results and a discussion. 



The constant p surfaces are coaxial ellipsoids. Projec- 
tion of these models on a plane perpendicular to a line of 
sight, and therefore, the calculation of ellipticity and po- 
sition angle are performed numerically. We refer to these 
models as M 2 models. 

Anoth er form of triaxial generalizat ion of (1) is inves- 
tigated by Ide Zeeuw and Carollol ( 1996), where two more 
terms are added to equation (1), each one of these is a suit- 
able radial function multiplied by spherical harmonics of low 
order. The models provide simple analytical representation 
of the observed surface brightness of triaxial elliptical galax- 
ies. However, for large values of flattening models become 
peanut shaped and are not used in the present investiga- 
tion. Very flat de Zeeuw - Carollo models are discussed in 
Appendix A. 



2 MODEL 



We use models, whic h are tr i axial generalizations of the 
spherical 7 models of iDehner] l| 19931 ). with density p given 
by 



p(r) = 



M (3- 7 )6 



4tt 



y (b + r) 



-4+7 



(1) 



where Mo is the mass of the model, r is the radial coordinate, 
^ 7 < 3 and b is the scale length. The models have cusp at 
the centre, and the density decreases as r~ at large radii. 
Dehnen's mode l s are the generalization of the well studied 
models of lJaffei (|l983T l and lHernquistl |l990), corresponding 
to 7 = 2 and 7=1, respectively. The projected surface den- 
sity of the model of Dehnen, corresponding to 7 = 1.5, most 
closely resembles to the de Vaucouleurs f? 1//4 law. Presently, 
we concentrate to 7 = 1.5 models only. 

A triax i al ge neralization of (1) is presented in 
IChakrabortvl (120041 ). which is modified in C08. The model 
is the density distribution of the same form as (1) with r 
replaced by M, where 



M 2 = x 2 



y 

pi 



with varying axial ratios 
p - 2 {M)= ^ 2 + M 2 p 



and 



\M) = 



13b 2 + M 2 
Pb 2 qo - 2 + M 2 q^ 2 



fib 2 + M 2 



(2) 



(3) 



(4) 



The axial ratios (P,Q) reduce to (po,<Zo) at small radii and 
to (j5oo,?oo) at large radii. /3 > is a parameter, which 
for a choice of (po, ?o,f>oo, goo) alters P and Q in the inter- 
mediate region. The models are fixed, once the axial ratios 
(po, 9o, Poo, goo) are chosen. The triaxialities To and Too are 
related to the axial ratios at small and at large radii by 



To = 



Pi 



T — 



1 



1 



(5) 



To fix up the scale length b of the triaxial models, we con- 
sider 7 = 1.5 and use the value of the effective radius 
R e = 1.286 of the spherical model. The effective radius of 
the triaxial models depends on the axial ratios, as well as 
on the vie wing angles. However, such changes are small for 
7 models |de Zeeuw and Carollol IT996). and are neglected. 



3 INTRINSIC SHAPES 

The galaxies chosen here are very flat. The morphological 
classification of NGC720, 2768 and 36 05 are E5/E3, E6/E5 
and EA/E5 respectively, from RC2 (|de Vaucouleurs et al.l 
Il976l ) catalogue. The apparent flattening of a elliptical 
galaxy depends on the intrinsic flattening and the orien- 
tation. Further, the marginal posterior density (MPD) V 
of the Bayesian estimate is obtained by integrating the pos- 
terior density over all viewing angles. To gain some insight 
into the possible values of the intrinsic shape, which will 
be obtained by Bayesian method, we perform the following 
numerical experiments. The objective of these experiments 
is to find suitable limits of the axial ratios (go, goo) for the 
plots of V . In the plots of V in C08, go and goo extend from 
0.5 to l.O. IStatlerl (|l994al ) chooses 0.4 < c L < 1.0. 

Fig. 1 shows the plot between the number N of viewing 
angles (9 ,(j> ) and the axial ratio q which gives ellipticity 
0.50 sC e < 0.55 (plot 1A) and 0.07 sC e sC 0.17 (plot IB). 
The number of viewing angles is counted between 0°.0 and 
90°. at the interval of 1°.0, both for 6 and tf> . The total 
number of viewing angles in this numerical experiment is 
8100. The axial ratio p is taken as 0.9. Here, we use Stark 
model. We find that a higher values of ellipticity is produced 
by flatter models and a lower values of ellipticity is produced 
by rounder models, over a larger number of viewing angles. 
Therefore, the Bayesian estimate should pick up a flat model 
to represent shape of the galaxies of our present investiga- 
tion. It is interesting to note that the plots (lA&zlB) show 
a maxima, which lies at q ~ 0.44 for the plot 1A and at 
q ~ 0.82 for the plot IB 

The fig. 1 and its inferences are based on applying Stark 
model, which has constant values of the axial ratio (p, g). 
However, in our shape estimates, we use models with varying 
axial ratios. So, we re-examine the results of these plots by 
considering M 2 models. 

Fig. 2 presents the marginal posterior density V as a 
function of (go, goo), summed over various values of (To, Too) 
for NGC720. We use the observational data, as shown in 
Table 1. We choose the values of both go and goo from the 
region 0.1 ^ (go,9oo) ^ 0.9. We use the M 2 models with 
P — 1.0. The probability of the shape is plotted in dark grey 
shade : darker is the shade, higher is the probability. The 
white contour encloses the region of 68% highest posterior 
density (HPD), which may be interpreted as la error bar. 
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Figure 1. Plot between the number of viewing angles and the 
axial ratio q, which would reproduce ellipticity in a chosen inter- 
val. Figure 1A is drawn for the ellipticity 0.50 ^ e ^ 0.55 while 
IB is drawn for the ellipticity 0.07 ^ e ^ 0.17. 



This figure indicates that higher probability region is con- 
fined, approximately between 0.3 to 0.8 of (go, Qoo)- Hence, it 
is more appropriate to choose the lower and the upper limits 
of both qo and qoo as 0.3 and 0.8 for the shape estimates of 
very fiat galaxies. This is discussed further, in sect. 3.1. 

For the same choice of the values of (qo,9oo), fig. 3 
shows shape - P(go,<?oo) of a rounder galaxy NGC3379. We 
use the observational data Ei„ = 0.078, e out = 0.133 at 
Rin = 15 .7 and R ou t = 49 .3. The effective radius R e of 
NGC3379 is 37" .5. We find that the HPD region is confined 
between(go , q<x) ^ 0.4 and the highest values of (50, 500) al- 
lowed in this plot. 



0.9 



ll.l 




0.1 ~ 0.9 



Figure 2. Plot of marginal posterior density (V) as a function 
of 90 , Qoo (= 9) , summed over various values of (To , Too ) , for NGC 
720 using the limits 0.1 to 0.9, both for and Qoo- 




Figure 3. Same as Fig 2, for NGC 3379. 



3.1 NGC720 

The observed data of NGC720, is taken from R-band sur- 
face photometry of IPeletier et all (|l99Qh . The ellipticity e 
increases monotonically from 0.315 at Ri„ = 8.5 arcsec to 
0.442 at Rout = 51.8 arcsec. In this range, the position angle 
decreases by 3°. 5. We consider the uncertainty in the ellip- 
ticity as 0.02 and in the position angle is 1°.0, both at Ri„ 
and at R ou t- These are the typical errors in observations 
(|de Carvalho et all Il99ll ; iPenereiro et all Il994h . The effec- 
tive radius of the galaxy is 52.0 arcsec. We use the ensemble 
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max=0.1997E-00 





0.5 



Figure 4. Plot of unweighted sum of MPD (P) as a function 
of qo,q<x(= q), for NGC 720 using the limits 0.5 to 1.0 both 
for go and <joo. The sum is taken over the M 2 models with = 
5.0, 2.5, 1.0, 0.5 and 0.2. Plus marks the location of the maximum 
probability. 





Figure 6. Three dimensional plot of the unweighted sum of MPD 
(V) as a function of qo, goo, \TJ\, for NGC 720. The sum is taken 
over the M 2 models with /3 = 5.0,2.5, 1.0,0.5 and 0.2. Values of 
\TJ\ arc constant in each section. In each section, qo goes from the 
left to right hand side from 0.25 to 0.75, and goo runs between 
the same values from the bottom to the top. 



0.80 



max=0.2523E-01 




0.30 



0.80 



Figure 5. Same as Fig 4, for NGC 720 using the limits 0.3 to 
0.8, both for qo and goo. 



of models, as described in sect. 2, with j3 = 5.0,2.5, 1.0,0.5 
and 0.2. Taking the sum of the marginal posterior density 
over all possible values of T and Too, and taking the un- 
weighted sum over all the models, we obtain shape estimate 
V as a function of (go, goo). 

Fig. 4 presents the shape estimate V (go, goo) of 
NGC720, wherein we have allowed the limits 0.5 to 1.0 both 
for go and goo- We find that the la region is very narrow 



which should be the consequence of the choice of the limits 
of go and goo. Examining this limit in fig. 1, we find that this 
choice falls in the region where high values of the ellipticity 
will not be reproduced. Therefore, we need to go to smaller 
values of (go, goo) to obtain higher ellipticities, which may 
be close to observed ellipticities of NGC720. 

Fig. 5 presents shape V (go, goo) wherein we have al- 
lowed the limits 0.3 and 0.8 for go and goo- If region is 
wider now (but narrow enough to satisfy the requirement 
of the likelihood dominated shape estimate). Although, it 
is the plot of MPD (V) as a function of shape parameters, 
which constitute the Bayesian estimate of the shape, some 
statistical summary of the shape is very convenient for its 
description. The expectation values < go >, < goo > and lo- 
cation of the peak values gop, gooP are such quantities. Table 
2 provides such a summary. The expectation values of the 
flattening at small and at large radii are < go >= 0.64 and 
< goo >= 0.43, respectively. 

Both in fig. 4 and 5, we choose the interval between 
higher and lower limits of (go, goo) as 0.5. This is basically 
to save the computer time, but maintaining the reliability 
of the results. Shape calculation requires a very large num- 
ber of projections, which need to be calculated numerically. 
We divide the parameter space of (go, goo) in 48 x 48 square 
bins of equal size and calculate the likelihood at the centre 
of each bin. The bin size is small enough, so that the calcu- 
lated likelihood can be regarded as a continuum function of 
(go, goo), and at the same time, the number of bins is small 
enough, so that the computer time is not unmanageable. 

Fig. 6 shows the 3-dimensional intrinsic shape of 
NGC720 as a function of go, goo and \T d \. We cut a total 
of 16 sections, each perpendicular to \T d \ axis, and arrange 
these sections in a form of a two-dimensional array. The 
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Table 1. Observational data of the galaxies. 



0,75 ■ 



max=0.2836E-01 



Galaxy 


R. e 




Rout 




tout 


e d 


NGC 720 


52.0 


8.5 


51.8 


0.315 


0.442 


-3.5 


NGC 2768 


76.5 


15.8 


95.4 


0.364 


0.569 


-1.4 


NGC 3605 


22.5 


5.9 


20.4 


0.305 


0.418 


-2.0 



Table 2. Statistical summary of the 2-dimensional shape esti- 
mates "P(go,goo) of the galaxies. 



Galaxy 


10p 


Qoop 


< 90 > 


< goo > 


AfGC720 


0.68 


0.48 


0.64 


0.43 


JVGG2768 


0.65 


0.29 


0.63 


0.32 


JVGG3605 


0.72 


0.49 


0.62 


0.42 



value of | | is constant in each section, and is shown in 
the plot. We find that the la region occupies larger area 
in the sections with smaller values of \Td\. Further, in each 
section of constant \Td\, la region occupies a small area of 
(9o 5 Qoo) plane. We find that higher V is concentrated in sec- 
tions with \Fj\ between 0.28 to 0.47. The expectation value 
of < IT* I >=0.41. 



3.2 Intrinsic shapes of NGC2768 and 3605 

The observational data used in the models of these galaxies 
are presented in Table 1. Here also, the data is obtaine d 
from R - band surface photometry of lPeletier et alj £1990) . 

Fig. 7 and 8 present the plot of V of NGC2768 and 
3605 as functions of (go, goo). The lower and upper limits of 
qo and goo are taken as 0.25 to 0.75. The HPD region shows 
that these galaxies are very fiat. The expectation values are 

< go >= 0.63 and < >= 0.32 for JVGG2768 and are 

< go >= 0.62 and < g^ >= 0.42 for 7VGG3605. We find 
that these galaxies are also intrinsically very fiat. Fig. 9 and 
10 present the 3-dimensional plot of V of NGC2768 and 
3605 as a function of go, goo and \T d \. 

Statistical summary of the intrinsic shapes of all the 
three fiat galaxies 7VGG720, 3605 and 2768 is presented in 
Table 3. Here, the values are taken from the 3-dimensional 
shape estimates. The expected and the peak values of go and 
goo as obtained from 2-dimensional estimates V (go, goo) are 
reported in Table 2. These values are quite close but not ex- 
actly the same as those reported in Table 3. The differences 
may be attributed to "resolution" : in 2-dimensional shape 
estimates, we have divided the parameters space of (go, goo) 
in 48 x 48 divisions, whereas in 3-dimensional shape esti- 
mate, space (go, goo) is divided in 10 x 10 divisions for each 
\T d \- 

Table 3. Statistical summary of the 3-dimensional shape esti- 
mates 7-" (go, goo, 1 2^ |) of the galaxies. 



Galaxy 


qo P 




\Td P \ 


< go > 


< goo > 


< \Td\ > 


NGC720 


0.68 


0.38 


0.41 


0.56 


0.40 


0.41 


NGC2768 


0.68 


0.28 


0.22 


0.62 


0.33 


0.39 


NGC3605 


0.68 


0.43 


0.16 


0.60 


0.41 


0.37 



0.25. 




0.75 



Figure 7. Same as Fig 4, for NGC 2768 using the limits 0.25 to 
0.75, both for go and goo- 



0.75 



max=0.1796E-01 




0.25 



0.75 



Figure 8. Same as Fig 4, for JVGG3605 using the limits 0.25 to 
0.75, both for go and goo- 



4 RESULTS AND DISCUSSION 

We have presented the intrinsic shapes of 3 very flat galax- 
ies. A specific feature of these estimates is the choice of the 
lower limit of go and g oo. The l ower li mit was chosen as 0.5 
by C08 and as 0.4 by IStatlerl (|l994al ibh in their investiga- 
tion of galaxies which are comparatively rounder, without 
mentioning any specific reason for this choice. Through the 
plots of fig. 1, and the detail discussion of the plots of V 
for 7VGG720, we justified the choice of very small values of 
(20, ?oo) as their lower limits for very flat elliptical galaxies. 
We took the lower limit as either 0.25 or 0.3. 
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Figure 9. Same as Fig 6, for JVGC2768. 
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Figure 10. Same as Fig 6, for NGCmOb. 



A summary of the intrinsic shapes of these very flat 
galaxies is presented in Table 2 and Table 3. We find that 
these galaxies are little rounder inside (average value of < 
qo >~ 0.6), but very flat outside (average value of < q^ >~ 
0.4). Following the nomenclature introduced in C08, these 
galaxies may be termed as RF type. 

Intrinsic shapes of elliptical galaxies have implications 
for their formation and evolution. As the galaxies studied 
here are very flat, we have given emphasis on the shape 
V{qo, qoo), and on informations about the flattening at inner 
and at outer radii. 



APPENDIX A: VERY FLAT DE ZEEUW - 
CAROLLO MODELS 

A simple family of triaxial models, with ellipticity 
variation and positi o n an gle twist was presented by 
Ide Zeeuw and Carollol fl996) with density distribution 

p(r, 9, <j>) = f(r) - g(r)Y 2 °(0) + h(r)Y ( f(9, 0) , (Al) 

where f(r) is same as (1), g(r) and h(r) are two radial func- 
tions, and y 2 ° and Y 2 2 are the usual spherical harmonics. 
Here, (r, 6, (j>) are the standard polar co-ordinates. The pro- 
jected surface density of (Al) can be calculated easily, and 
often analytically. g(r) and h(r) are fixed by assiging ax- 
ial ratios (po,9o) and (poo,9oo) respectively, at small large 
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Table Al 


. Regions of negative 


P- 


P = PO = 


Poo =0.9,7 = 1.5,9 


= o°.o,4>= f. 


q = qo = q. 


now 


r high 
b 


0.55 


1.93 


> 6.40 


0.56 


2.31 


6.27 


0.57 


3.08 


4.48 


^ 0.58 


p is positive at all r 



radii, where constant - p surfaces are approximately ellip- 
soidal. Numerical distribution function was shown to exist 
for prolate triaxials : (p, q) = (0.65, 0.60) and for oblate 
triaxials : (p, q) = (0.95,0.65), wherein it is assumed that 
go = goo = g,po — p<x = p and 7 = 1.0 or 1.5 (jThakur et al.l 
l2007h . The rounder versions of these models are employed 
successf ully in many investigations, incl uding the shape es- 
timates (|Thakur and ChakrabortvfeoOll ; C08). 

However, very flat versions of de Zeeuw - Carollo mod- 
els have several unde s irable features. It was realized by 
Ide Zeeuw and Carollol dl996h that the constant p surfaces 
become peanut shaped or dimpled for large values of flat- 
tening. Such models can not be a "true" representation of 
the shape of an elliptical galaxy. 

In addition to above, we now find the appearance of 
narrow regions, where p is negative. Clearly, it is unphysical 
to call such negative p as mass density. Such negative p 
appears in polar regions (6 ~ 0°.0), at an intermediate r 
extending from some ri ow to rhigh- Tables (Al) and (A2) 
show some of the regions of negative p on cj> = ? plane. 

Triaxials models of similar form as (Al) was proposed 
by ISchwarz schildl (|l979l ) as a numerical model, where f(r) 
is taken as the modified Hubble dens ity distribution. Later, 
it was put into an analytical form by Ide Zeeuw and Merrittl 
(1983). Projected propert ies of such triaxial modified Hub - 
ble model was studied bv lChakrabortv and Thakurl (|2000h . 
We now find that sufficiently flat versions of triaxial modi- 
fied Hubble model also exhibit regions of negative p. 

We find that the appearance of negative p regions is 
correlated with the dimpleness of constant p surfaces. We 
examine an extended version of models (Al), which includes 
terms with high order spheri cal harmonics Y® ' ,Y% and Y4 . 
Such models were studied bv lChakrabortv and Dai (|2003h . 
It was found that the dimpleness reduces i.e., the models 
become more ellipsoidal - like. We now find that for the 
same choice of the parameters as in Table (Al), the interval 
(rhigh - n ow ) of negative p decreases. 

It will be interesting to extend the studies of particle or- 
bit an d the numerical distribution function of lThakur et al.l 
(120071 1. to the models which exhibit regions of negative p. 

This paper has been typeset from a TfrjX/ WFftX. file prepared 
by the author. 



Table A2. Regions of negative p. 



P = P0 


= Pao =0.9,7 = l-5,</>= § ,g = 0.57 
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o°.o 


3.08 


4.48 


1°.0 


3.08 


4.48 


2°.0 


3.21 


4.36 


3°.0 


3.46 


3.85 


4°.0 


p is positive at all 
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